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ABSTRACT 

Far-infrared Spitzer observations of elliptical galaxies are inconsistent with simple steady state 
models of dust creation in red giant stars and destruction by grain sputtering in the hot interstellar 
gas at T ~ 10 7 K. The flux at 24/im correlates with optical fluxes, suggesting that this relatively 
hot dust is largely circumstellar. But fluxes at 70/zm and 160/xm do not correlate with optical fluxes. 
Elliptical galaxies with similar Lb have luminosities at 70/zm and 160/um (X70 and £160) that vary 
over a factor ^100, implying an additional source of dust unrelated to that produced by ongoing local 
stellar mass loss. Neither Ljq/Lb nor Liqq/Lb correlate with the stellar age or metallicity. Optical 
line fluxes from warm gas at T ~ 10 4 K correlate weakly with L70 and Liqq, suggesting that the dust 
may be responsible for cooling this gas. Many normal elliptical galaxies have emission at 70/im that is 
extended to 5-10 kpc. Extended far-infrared emission with sputtering lifetimes of ~ 10 s yrs is difficult 
to maintain by mergers with gas-rich galaxies. Instead, we propose that this cold dust is buoyantly 
transported from reservoirs of dust in the galactic cores which is supplied by mass loss from stars in 
the core. Intermittent energy outbursts from AGNs can drive the buoyant outflow. 
Subject headings: galaxies: elliptical and lenticular; galaxies: ISM; infrared: galaxies; infrared: ISM 



1. INTRODUCTION 

The study of far-infrared (FIR) emission from ellipti- 
cal galaxies is particularly interesting because the inter- 
nal rates that dust is formed and destroyed can both be 
estimated. Dust is produced by mass loss from evolv- 
ing red giant stars in the old stellar population typical 
of these galaxies. This mass loss is directly observed 
in the cumulative mid- infrared 10-24/mi emission from 
circumstellar regions around giant stars, causing the sur- 
face brightness profiles at 10-24/mi to resemble those at 
optical wavelengths (Knapp et al. 1989; Athey et al. 
2002). The dust destruction rate by sputtering can be 
estimated from the frequency of thermal impacts from 
ions in the hot interstellar gas which for many ellipti- 
cals can be estimated from X-ray observations. During 
their sputtering lifetimes, dust grains ejected from gi- 
ant stars are heated by starlight and thermal electron 
impacts to temperatures that radiate in the FIR (e.g. 
Tsai & Mathews 1995, 1996; Temi et al. 2003). This 
radiation has been observed both by the Infrared Astro- 
nomical Satellite, IRAS (Knapp et al. 1989) and by the 
Infrared Space Observatory, ISO (Temi et al. 2004). In 
addition, small optically obscuring clouds and disks of 
dusty gas have been detected in about half of the bright 
ellipticals studied with the Hubble Space Telescope, HST 
(van Dokkum & Franx 1995; Ferrari et al. 1999; Lauer 
et al. 2005). These clouds probably formed from dusty 
gas ejected from stars within the central ~kpc in these 
elliptical galaxies (Mathews & Brighenti 2003). 
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Galaxy mergers are another obvious possible source of 
dust in elliptical galaxies, indeed it is widely thought 
that ellipticals are formed by mergers. Elliptical- looking 
galaxies can be produced after a few Gyrs following the 
merger of two disk galaxies (e.g. Barnes & Hernquist 
1992), but this simple recipe may not easily explain 
all structural and kinematic properties of these galax- 
ies (Naab, Jesseit & Burkcrt 2006). Many of the merg- 
ers that produced the ripples, shells and counter-rotating 
cores in elliptical galaxies probably occurred a very long 
time ago. It is likely that massive boxy ellipticals, sim- 
ilar to many of those with strong FIR emission, have 
experienced mostly "dry" gas and dust-free mergers dur- 
ing the last few Gyrs (Bell et al. 2004; Faber et al. 
2006). Nevertheless, mergers do occasionally occur be- 
tween low-redshift ellipticals and gas-rich galaxies and 
these post-merger ellipticals often selectively appear in 
the observing lists of infrared observers. Our previous 
study of FIR emission from elliptical galaxies was based 
on a sample chosen from the ISO data archive (Temi et 
al. 2004) that included a significant number of galaxies 
with abnormally large masses of dust and cold gas which 
are uncharacteristic of elliptical galaxies in general. The 
large range of FIR luminosities that we found in the ISO 
sample was therefore not unexpected, and it was rea- 
sonable to suppose that mergers with dusty galaxies ex- 
plained the absence of correlations between optical and 
FIR emission. The failure to find correlations between 
the FIR and optical luminosities of elliptical galaxies has 
been used in the past to argue that the cold dust is ac- 
quired in a stochastic fashion by mergers (Forbes 1991). 

To understand the relative role of internal and external 
sources of dust in elliptical galaxies, we describe here new 
Spitzer observations at 24, 70 and 160/im. Our full sam- 
ple consists of 46 elliptical galaxies, formed by combin- 
ing our own proprietary Spitzer observations with those 
of other observers. The improved sensitivity and angu- 
lar resolution of Spitzer provides clear evidence that the 
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70/im emission is spatially extended in many ellipticals. 
In particular, we find extended FIR emission in ellipti- 
cal galaxies that do not contain uncharacteristically large 
masses of neutral, molecular or optically visible cold gas. 
The diffuse nature of the dust in these otherwise normal 
galaxies, combined with the relatively short (ion sputter- 
ing) lifetime for dust destruction in the hot gas, favors 
an internal source of dust rather than dust-rich merg- 
ers which are probably too infrequent. The luminosity 
of extended FIR dust emission in most ellipticals greatly 
exceeds that expected from normal mass loss from local 
stars. However, small, optically thick nuclear dust clouds 
observed in many elliptical galaxies may provide an alter- 
native, more plausible source for the extended FIR dust 
emission. Consequently, we argue that dust accumulated 
in the central regions of elliptical galaxies from stellar 
sources in the core is transported out to larger radii 
by buoyantly rising gas intermittently heated by low- 
luminosity active galactic nuclei (AGN) in the galactic 
cores. If this hypothesis can be confirmed in more detail, 
the presence of extended cold dust can be used to trace 
the radial circulation of hot gas in these galaxy /groups. 
Such circulation flows are known to be consistent with 
X-ray observations in which very little hot gas cools to 
low temperatures (Mathews, Brighenti & Buote 2004). 

2. THE SAMPLE 

For this study we construct a sample of 46 elliptical 
galaxies by combining far-infrared Spitzer observations 
from our own Spitzer data with additional galaxies in 
the public Spitzer archive. Our proprietary far-infrared 
Spitzer observations of elliptical galaxies (PID 20171, PI 
P. Temi) consists of 11 non-peculiar, bright (Lb > 2xl0 10 
Lb,q) elliptical galaxies, most having high-quality X-ray 
observations. Galactic proximity and spatial scale were 
also important factors in choosing this sample since high 
angular resolution is needed to separate emission from 
dusty cores and diffuse emission from truly interstellar 
dust spread over the galactic volume. Furthermore, we 
avoided sources with strong radio sources and AGN that 
have detectable signals in the far infrared. Ellipticals 
with strong infrared Milky Way cirrus or bright spiral 
galaxies in the vicinity were also avoided. Since one of 
our goals is to study the interaction between dust and 
hot gas, we included many galaxies for which the X-ray 
emission is well-studied by either Chandra or XMM. 

In addition to our original sample we selected 35 ad- 
ditional ellipticals from data currently available in the 
public Spitzer archive targeted by several Pis during their 
General Observation in Cycles 1 and 2, Guaranteed Time 
Observations, and the Legacy Program SINGS (Kenni- 
cutt et al. 2003). All the galaxies added to our original 
sample were chosen using selection criteria tailored to 
our specific scientific goals. Thus, the 46 galaxies stud- 
ied here form a rather homogeneous sample and is large 
enough to be statistically meaningful. 

Tables 1 and 2 summarize the relevant properties of 
the sample galaxies including the far-infrared photomet- 
ric measurements. Optical data are taken from the 
LEDA catalog (Paturel et al. 1997) and the RSA catalog 
(Sandage & Tamman 1981), while X-ray data are from 
Ellis and O 'Sullivan et al. (2006). Luminosities and dis- 
tances are calculated with H = 70 km s _1 Mpc -1 and 
are derived from the flux densities in Table 1 using far- 



infrared (FIR) passbands of AA = 5.3, 19 and 34/im for 
emission at 24, 70 and 160/im respectively. The Pis of 
the original Spitzer program from which the galaxies have 
been selected are shown in column two of Table 1. All 
galaxies are classified as ellipticals in the NED catalog, 
and the morphological type T is taken from the LEDA 
catalog. Nevertheless, we shall find that our sample con- 
tains some SO galaxies. Most of the sample galaxies have 
known X-ray luminosities L x , which are a measure of the 
hot gas present in each galaxy. 

Of the 46 galaxies in the sample three are dwarf ellip- 
tical galaxies (NGC855, NGC1377, and NGC3265) with 
Lb < 5 x 10 9 Ls !0 . Although they are listed in Table 1, 
they are not included in our discussion. 

3. OBSERVATIONS AND DATA REDUCTION 

The data presented here were acquired with the Multi- 
band Imager Photometer (MIPS) (Ricke et al. 2004) 
on the Spitzer Space Telescope (Werner et al. 2004) in 
three wavebands centered at 24, 70 and 160/xm. Since the 
galaxy sample has been constructed by collecting archival 
data from several Spitzer observing programs, the image 
data were recorded using observing modes and integra- 
tions times tailored for their original scientific goals. As 
a consequence, although the maps represent the deepest 
FIR images acquired to date, some data achieve sensitiv- 
ities that are not always optimized for targets that are 
intrinsically faint at far-infrared wavelengths. For exam- 
ple, three of the brightest ellipticals in the Virgo Cluster 
(NGC4472, NGC4636, and NGC4649) recorded under a 
guaranteed time program (PID 69, PI G. Fazio) reach 
a relatively low sensitivity of 0.5MJy/sr and l.lMJy/sr 
(lcr) at 70 and 160^to, respectively, while other galax- 
ies (i.e. PID 20171, PI P. Temi) have deeper maps at 
a sensitivity level of only 0.12MJy/sr and 0.3MJy/sr for 
the same two wavebands. The SINGS data are recorded 
in MIPS scan mode, covering a very large sky area 
(~ 30' x 10'), incorporating two separated passes at each 
source location. SINGS images correspond to maps with 
intermediate sensitivity. Details on the observing strate- 
gies, field coverage and integration times for the SINGS 
program are described by Kennicutt et al. (2003). Apart 
from the SINGS observations, data have been acquired 
in MIPS photometry mode, allowing appropriate cover- 
age of the sources and their extended emission. At all 
wavelengths, sufficiently large field sizes were chosen to 
reliably measure the background emission. We used the 
Basic Calibrated Data (BCD) products from the Spitzer 
Science pipeline (version 13.2) to construct mosaic im- 
ages for all objects. Pipeline reduction and post-BCD 
processing using the MOPEX software package (Makovoz 
et al. 2006) provide all necessary steps to process indi- 
vidual frames: dark subtraction, flat-fielding, mux-bleed 
correction, flux calibration, correction of focal plane ge- 
ometrical distortion, and cosmic ray rejection. 

Foreground stars and background galaxies were present 
in the final mosaiced images at all bands. These were 
identified by eye and cross-checked using surveys at other 
wavelengths (Digital Sky Survey and 2MASS). They were 
then masked out from each MIPS image before flux ex- 
traction was performed. Also, a few asteroids were iden- 
tified and removed from the affected maps. Aperture 
photometry was performed by measuring fluxes in suit- 
able circular apertures around the centers. Flux densities 
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were extracted from apertures extending out to the entire 
optical disk of each galaxy (i?2s)- Sky subtraction was 
performed by averaging values from multiple apertures 
placed around the target, avoiding any overlap with the 
faint extended emission from the galaxy. Uncertainties 
in flux densities due to sky subtraction are generally less 
than 1%, but can be of the order of tens of percent for 
faint sources. 

The measured flux densities are listed in Table 1. Most 
of the galaxies in the sample have been detected by 
MIPS, and among them a large fraction show far-infrared 
emission well beyond the sensitivity capabilities of pre- 
vious infrared space observatories. Although most of the 
galaxies presented in this paper are unpublished either 
because they are our own proprietary data or because the 
analysis is still in progress by the original Pis, results 
for a few galaxies are already present in the literature 
(i.e. the SINGS galaxies). We find that our flux den- 
sities are in good agreement with previously published 
data. The only exception is the measured flux density 
at 160^m reported by Dale et al. (2005) for the SINGS 
galaxy NGC584. Very likely this discrepancy is due to a 
misidentification of the target in their very large 160^m 
scan map. Corrections for extended emission were ap- 
plied to the fluxes as described in the Spitzer Observer's 
manual. Photometric errors listed in Table 1 refer only to 
the statistical uncertainties. Total errors should include 
systematic uncertainties due to flux calibration and to 
sky subtraction errors. The uncertainties on the final 
absolute calibration are estimated at 10% for the 24/xm 
data, and 20% at 70 and 160^m. 

Figure 1 shows gray-scale images of a representative 
sub-sample of galaxies. The column of panels on the left 
is the 24^m data, while 70 and 160^im data are shown in 
the center and right columns, respectively. Galaxies are 
ordered by ascending IR flux in a common field of view 
of 3' x 3'. Some very faint galaxies are shown, including 
a non detection at 160/xm and few sources with extended 
far-infrared emission. Detailed images of the diffuse ex- 
tended emission in NGC4636 and NGC5044 will be pre- 
sented in a companion paper currently in preparation. 

4. COMPARISON OF FIR, OPTICAL AND X-RAY EMISSION 

In Figure 2 we compare fluxes at 24, 70 and 160^m 
with optical B-band fluxes and a related plot of FIR and 
B-band luminosities is shown in Figure 3. Both the 24/xm 
and B-band fluxes and luminosities are clearly correlated, 

logF 24 = (8.37 ± 0.42) + (0.716 ± 0.043) logF B , Q (1) 

and 

\ogL 2A = (32.37 ±0.69) + (0.799 ±0.066) log L B;0 . (2) 

The approximate consistency with unit slope F 2A oc Fb 
indicates that the emission at 24/im, like that at mid-IR 
wavelengths (Knapp et al. 1989; Athey et al. 2002), is 
dominated by circumstellar emission that closely traces 
optical starlight. 

However, the huge scatter seen in the 70/xm-Bband 
panels in Figures 2 and 3 suggest that this colder dust 
is unassociated with stars and is truly interstellar. The 
absence of any clear correlation between FIR and opti- 
cal emission from ellipticals has been discussed by Forbes 
(1991), Goudfrooij & de Jong (1995), Trinchieri & Goud- 
frooij (2002) and Temi et al. (2004). In principle cold, 



FIR-emitting dust could either be diffusely distributed 
in the hot gas or confined to dense central clouds similar 
to those visible in optical absorption (e.g. Lauer et al. 
2005). We discuss below our preference for the diffuse 
component in producing most of the FIR-cmission. 

It is remarkable that L70 and L160 range over more 
than two orders of magnitude for galaxies with similar 
Lb - this is a significant new result that must be ex- 
plained. In our previous survey of archival data from In- 
frared Space Observatory (ISO) we found a similar spread 
in FIR luminosities at 60-200/im, but the FIR luminosi- 
ties in the ISO data spanned only one order of magni- 
tude for fixed Lb (Temi et al. 2004). At the time we 
attributed much of this scatter to the unusual nature of 
the ISO sample which contained many (post-merger) el- 
liptical galaxies with massive disks of cold, dusty gas. 
Now with this Spitzer sample of (mostly) normal ellipti- 
cal galaxies we find that the range of FIR luminosities is 
even larger, although we show in Appendix A that some 
of the FIR-brightest sample galaxies cannot be regarded 
as normal. In Figure 4 we show that there is no obvious 
trend of dust color F^q/Fiqq (i.e. mean dust tempera- 
ture) with L70 except perhaps a tendency for cooler dust 
in more FIR-luminous galaxies. A similar conclusion re- 
sults if we replace F70 on the horizontal axis with i<i6o, 
L70 or Li6o- 

Figures 5 and 6 show no discernable correlation of ei- 
ther L^o/Lb or Liqo/Lb with stellar Balmcr line ages 
or metallicity. However, L 2A /Lb may show a slight ten- 
dency to decrease with age 

logL 2i /L B = (30.35 ± 0.05) - (0.021 ± 0.005)age (3) 

or increase with metallicity 

logL 2i /L B = (29.95 ±0.06) + (0.7 ± 0.2) log[Z/H\. (4) 

A negative correlation of L 2 ^/Lb with age would be ex- 
pected since the rate of mass loss from evolving stars 
decreases with stellar age (Piovan et al. 2003), although 
this expected trend is not supported by our study of ellip- 
tical galaxies at mid-IR wavelengths (Temi et al. 2005). 
It must be noted however that Balmer line ages for the 
same galaxy often show considerable variation when de- 
termined by different authors, so the anticorrelation of 
L 2 a/Lb with age needs further confirmation. A posi- 
tive correlation of L 2 ±/Lb with metallicity in Figure 6 
is expected if more interstellar grains are produced by 
mass-losing stars of higher metallicity. However, this is 
inconsistent with the negative correlation at 70/xm in Fig- 
ure 6, 

log L 70 / L B = (30.7 + 0.1) - (1.9 ± 0.4) log[Z/H]. (5) 

No correlations with X-ray luminosity with any of the 
FIR luminosities are visible in Figure 7. Decreasing FIR 
luminosities with increasing hot gas density might be ex- 
pected since interstellar grains are destroyed by sputter- 
ing which is proportional to the density of the hot gas. 
However, the relationship between L x and gas density 
is unclear. For example, there is no correlation between 
L x and the hot gas density at ~ 10 kpc in the ellipti- 
cal galaxies observed with Chandra by Humphrey et al. 
(2006). This disconnect between gas density and L x may 
arise because the hot gas density on kpc-galactic scales 
is influenced by intermittent energy supplied by the cen- 
tral black holes, while L x is largely determined by more 
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extensive gas that fills the dark halos that surround the 
elliptical galaxies (Mathews et al. 2006). Consequently, 
the absence of correlations in Figure 7 is not surprising. 
Finally, we note that all qualitative features in Figure 7 
are preserved if luminosities are replaced with fluxes. 

In Figure 8 we compare FIR fluxes with optical emis- 
sion line fluxes determined by three different observers: 
Goudfrooij et al. (1994a,b), Macchetto et al. (1996) and 
Sarzi et al. (2006). The optical observations require dif- 
ficult background subtraction and the procedures used 
vary among the observers. For this reason, we compare 
FIR fluxes in Figure 8 with each set of optical emission 
line fluxes independently 

We believe that both F 60 and F170 correlate weakly 
with observed optical line fluxes, but that F24 does not. 
The recent SAURON data for F Hf3 (Sarzi et al. 2006) 
show the tightest relationship. The Goudfrooij et al. 
data (left column of panels in Figure 8) are marginally 
correlated, but the significance increases if it is combined 
with the Macchetto et al. data (shown in the central col- 
umn). 

In Figure 9 we compare the fluxes in the optical emis- 
sion lines with B-band fluxes, using only data for non- 
peculiar E or E/S0 galaxies, excluding peculiar and 
known SO galaxies. The top and bottom panels from 
Macchetto et al. (1996) and Sarzi et al. (2006) show 
that the line flux is weakly correlated with the stellar 
continuum (see also Phillips et al. 1986; Macchetto et al. 
1996; Goudfrooij 1999). A correlation between extended 
line emission and the optical continuum suggests that 
the warm gas with T ~ 104 K has a stellar origin. But 
the quality of these difficult optical observations is poor. 
There are very few objects in common among the Goud- 
frooij, Macchetto and Sarzi data sets, but the flux- flux 
plots of these common galaxies show large scatter with 
almost no clear monotonic relationship. Five of the six 
galaxies in common between Goudfrooij and Macchetto 
correlate negatively in a flux- flux plot. 

The correlation in Figure 9, if real, may simply indicate 
that the warm gas and the dust associated with it both 
have a stellar origin - more stars of a similar age and 
mctallicity eject more gas. However, the optical emis- 
sion line flux depends sensitively on the density of the 
warm gas and therefore on the pressure in the local hot 
interstellar gas. At present there are too few accurate 
observations to interpret Figure 9 with confidence. 

All ellipticals in our sample show evidence of extended 
emission at 24/mi similar to that of the optical stellar 
images - details of this will be discussed in a separate 
publication. However, spatial extension is also evident 
at 70^m as shown in Table 1, but the emission at 160^m 
is not spatially resolved due to decreasing resolution at 
larger wavelengths. Figure 10 shows the surface bright- 
ness at 70/^m for eight galaxies together with the point 
response function (PRF) for the MIPS detector (includ- 
ing the first Airy ring). All contours in Figure 10 are 
arbitrarily normalized to 50 at the center. The (loga- 
rithmic) extension in arseconds beyond the PRF is best 
measured at 0.08 of the central fluxes, shown with arrows 
in Figure 10. The full width at 4/50 = 0.08 maximum 
(subsequently referred to as FW0.1M; last two columns 
in Table 1) of the PRF is 37". Twelve galaxies in our 
full sample have FW0.1M greater than 40", indicating 
that these galaxies are spatially resolved at 70^m. When 



FW0.1M/2i? e is plotted against F^/Fi 60 , there is no 
clear evidence that galaxies with more extended F70 have 
colder dust. 

5. A SIMPLE FAR-INFRARED SPECTRAL ENERGY 
DISTRIBUTION 

We estimate the FIR spectral energy distribution 
(SED) from diffusely distributed interstellar dust in el- 
liptical galaxies using a simple model similar to that de- 
scribed in Temi et al. (2003). For completeness and 
clarity, we briefly describe this model again here. We as- 
sume that dust grains ejected from evolving stars come 
into direct contact with the hot interstellar gas after a 
short time, < 10 5 — 10 6 yr, when the stellar gas has 
heated to the (virial) temperature of the hot interstellar 
gas. This time scale is established by comparing typical 
Balmer emission line observations with the luminosity 
expected if all the gas expelled from stars (at the ex- 
pected rate for an old stellar population) is ionized and 
in hydrostatic equilibrium with the hot gas (Mathews & 
Brighenti 1999). After ~ 10 6 yrs this warm gas at 10 4 K 
must merge with the hot gas since otherwise its Balmer 
line luminosity would exceed typically observed values. 
However, the dust in this gas can resist sputtering de- 
struction in the hot interstellar gas (T ~ 10 6 — 10 7 K) for 
much longer times, ~ 10 8 yrs. Consequently, we model 
"recently" produced dust that was created by local stars 
during the previous ~ 10 8 yrs. Once in contact with the 
hot gas, the dust grains are destroyed by sputtering at a 
rate 
da 
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where a is the grain radius (in /zm) and n p is the proton 
density in the hot gas (Draine & Salpeter 1979; Tsai & 
Mathews 1995). 

If N(r,a,t)da be the number of grains cm' 3 with ra- 
dius between a and a + da at galactic radius r at time t, 
the grain population evolves according to 
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where S(r, a)da is the rate at which grains with radii a 
are expelled from stars per cm 3 (Tsai & Mathews 1995). 
We assume S(a) = S a~ s for < a < a max , with s = 3.5 
(Mathis et al. 1977), so the coefficient 
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depends on the rate that the old stellar population ejects 
mass a* = M~ 1 \dM*/dt\ « 4.7 x 10~ 20 s' 1 (Mathews, 
1989), the stellar density p*, the initial dust to gas mass 
ratio scaled to the stellar metal abundance in solar units 
5 = (l/150)z(r), and the density of (amorphous) silicate 
grains, p g = 3.3 gm cm~ 3 . The solution of eq. (7) is 

-1 
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The grain temperature Td is determined by the balance 
of heating (by both absorption of stellar radiation and by 
electron-grain collisions) and radiative cooling: 
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= 4na 2 c7sBT d (r, a) 4 {Q abs )(T d , a). (10) 

In view of the high electron-grain collision frequency, we 
neglect stochastic heating due to single impacts. 

We assume that all grains have properties similar 
to astronomical silicates (Laor & Draine 1993) using 
the following approximations for the absorption coeffi- 
cient Q a b s ~ atp(X), and for the Planck-averaged values 
(Qabs) ~ 1-35 x 10" 5 Tja (see Draine & Lee 1984). The 
mean intensity of starlight is found by integrating over 
a de Vaucouleurs stellar profile using the appropriate ef- 
fective radius R e and stellar luminosity for each galaxy. 
The collisional heating term in the equation above in- 
cludes a correction r(a) for small grains that do not com- 
pletely stop the incident thermal electrons (Dwek 1986) 
and (v e E e ) = (32/7rm e ) 1 / 2 (fcT) 3 / 2 is an average over a 
Maxwellian distribution. The electron density and tem- 
perature of the hot gas are calculated using dcconvoluted 
profiles from X-ray observations. The FIR emissivity at 
each wavelength is calculated from 

j(r,\) = ±- f N(r,a)4ira 2 l0- s Q abs (a,\) 
4?r Jo 

x TrB(T d (r,a), X)da. (11) 

In Figure 11 we compare the expected FIR spectral 
energy distribution for three galaxies in our sample that 
have well-known hot gas density and temperature pro- 
files from X-ray observations: NGC 1399, 1404 and 4472. 
None of these galaxies are spatially resolved at 70^m. 
The X-ray data for these galaxies has been taken from the 
following sources: NGC 1404 & 1399: (Paolillo et al 2002; 
Jones et al. 1997); NGC 4472: (Mathews & Brighenti 
2003; Brighenti & Mathews 1997). The profiles in Fig- 
ure 11 show the expected FIR spectral energy distribu- 
tions computed with the equations above. The solid and 
dashed lines correspond respectively to emission from 
grains with original maximum sizes a max = 0.3/im and 
l.O^m. The observed fluxes for NGC 1399 and NGC 4472 
at 70/xm and 160/xm are in excellent agreement with our 
simple model. The fit to NGC 1404 is also satisfactory, 
given the many assumptions in our model. While the 
flux at 24/xm is also shown in Figure 11, this circumstcl- 
lar emission is not included in our model, but it supports 
our assumption that much of the diffuse interstellar dust 
in these three galaxies has a stellar origin. 

It is of particular interest to determine FIR spectral 
energy distributions for galaxies that have accurate in- 
terstellar gas density and temperature profiles and are 
also spatially resolved at 70/xm. Seven galaxies in Ta- 
ble 1 have FW0.1M greater than 50" but only six of 
these are in the Ellis and O'Sullivan (2006) list of X-ray 
galaxies detected with ROSAT: NGC 2768, 2974, 4125, 
4636, 4696, and 5044. Of these six galaxies, only three 
- NGC 4636, 4696 and 5044 - have X-ray luminosities 
(log 2^ > 41.5) large enough for an accurate determina- 
tion of the gas density and temperature profiles. How- 
ever, the optically peculiar elliptical NGC 4696 is the 
central galaxy in the Ccntaurus Cluster with powerful 
FRI-type radio emission, a large mass of dusty neutral 
gas and highly disturbed hot cluster gas (Fabian et al. 
2006). In view of these complexities, NGC 4696 is an 
inappropriate galaxy to learn about the nature of FIR 
emission from normal elliptical galaxies. We find that 



the observed fluxes for the two remaining FIR-extended 
galaxies - NGC 5044 and NGC 4636 - are very much 
larger than predicted by our steady state model for the 
FIR spectral energy distribution. We believe that the 
dust in these galaxies has an entirely different origin. 
The large and extended FIR emission from NGC 5044 
and NGC 4636, which is related to the enormous range 
of FIR luminosities in Figure 3, will be discussed in a 
separate forthcoming publication. 

6. DISCUSSION 

6.1. Large Range of LpjR 

The first important conclusion to be drawn from Fig- 
ures 2 and 3 is that warm circumstellar dust and colder 
interstellar dust behave very differently. The relatively 
small scatter in F24 in Figure 2 is slightly larger than our 
estimated observational errors. Nevertheless, we have 
been unable to find correlations between stellar age or 
metallicity and the residuals in F24 from the mean cor- 
relation. We conclude that the warm dust emitting at 
24^tm is largely circumstellar and this is confirmed by the 
similarity of optical and 24/im surface brightness profiles 
which will be discussed in detail elsewhere. 

The second and most remarkable result in Figures 2 
and 3 is that the emission from colder dust at 70 and 
160/im varies by a huge factor relative to that at 24/im 
and optical wavelengths. In the previous section and 
in Figure 11 we show that 70 and 160/xm fluxes from 
three galaxies having well-observed interstellar gas den- 
sity profiles can be fit reasonably well with our simple 
dust heating model without varying any of the dust pa- 
rameters from those we assumed at the outset. 

However, notice that all three galaxies in Figure 11 
that agree with our model SED - NGC 1399, 1404 and 
4472 - have values of L70 and Li 60 in Table 2 that lie 
near the bottom of the observed distributions in Figure 
3. We think it is impossible to explain the huge disper- 
sion of L70 and Liqq in Figure 3 by simply varying the 
dust parameters in the SED models used for the three 
galaxies in Figure 11, i.e., the stellar mass loss rate, dust 
composition, dust size distribution, etc. are unlikely to 
be in error by an order of magnitude or more. Our SED 
model clearly fails to fit most of the FIR observations. 

Some of the large spread of L70 and Li &0 in Figure 
3 can be understood because our sample includes very 
FIR-luminous elliptical galaxies that contain abnormal 
masses of interstellar dust and cold gas, probably ac- 
quired by a recent merger. In Appendix A we review 
relevant observations of 12 sample galaxies having the 
highest 70^iin luminosities in order of descending L70. 
We conclude from this limited information that most el- 
lipticals with the largest L70 and L160 in Figure 3 cannot 
be regarded as typical or normal. Of the 12 galaxies in 
Appendix A, we regard the following 8 as either likely SOs 
(which often have extended, rotationally supported, star- 
forming disks of cold dusty gas) or having massive HI or 
molecular lanes or disks probably acquired by mergers: 
IC 3370, NGC 5018, 2974, 4125, 5322, 2768, 3962, and 
5077. These atypical or mis-classified elliptical galaxies 
are the brightest galaxies in our sample. However, we re- 
gard the following galaxies in Appendix A with slightly 
lower L70 as normal or quasi-normal ellipticals: NGC 
3557, 5044, 4589, and (possibly) 4697. If we disregard 
the eight galaxies with atypical dust content, the range 



6 



Temi, Brighenti & Mathews 



in L70 and L\qq occupied by normal ellipticals is reduced 
from ~ 100 (shown in Figure 3) to about ~ 30, but this 
is still larger than can be explained by varying dust pa- 
rameters in our SED model. These normal ellipticals 
with excess FIR are the most interesting galaxies in our 
sample. 

6.2. Location of FIR- Emitting Dust 

We assume that galaxies with Lfir significantly above 
that predicted by our model SEDs in Figure 11 have 
dust contents that exceed those that can be understood 
by a steady state balance between dust creation by mass- 
losing stars and dust destruction by sputtering. Galaxies 
in Figure 3 with logi7o <; 40.5 must contain additional 
cold dust from another source. This additional cold dust 
is either (1) spread diffusely throughout a large region of 
the galactic hot interstellar gas or (2) located in dense 
optically thick clouds similar to those seen in absorption 
by optical observers. In Appendix B we describe a sim- 
ple argument that supports the first option above, diffuse 
FIR emission. We suppose in Appendix B that the FIR- 
emitting dust is in a small spherical cloud of radius r c 
located at the center of an elliptical galaxy. The dust is 
heated both by internal stars and by starlight incident 
from the entire galaxy. The dust cloud is assumed to be 
optically thick to most stellar radiation. As the cloud ra- 
dius r c is increased, the estimated 70^m luminosity also 
increases because the cloud encloses more stars and in- 
tercepts more light from external stars. However, this 
simple dusty cloud model cannot account for the lumi- 
nosity range in Figure 3 unless r c is so large (> 200 pc) 
that it would be easily apparent in moderate resolution 
optical images. The UV-ionized surfaces of these clouds 
may also be visible. While AGN heating within optically 
thick dusty clouds could easily provide the needed FIR 
energy, this seems unlikely since AGNs would need to 
be bolometrically brighter than the stars within r c and 
would therefore be seen in those ellipticals in which the 
dusty clouds are in disks or not present at all. Such lu- 
minous AGN would almost certainly be visible at other 
frequencies. 

Comparisons of the most FIR-luminous (normal) ellip- 
ticals in our sample with recent HST images (e.g. van 
Dokkum & Franx 1995; Lauer et al. 2005), show a corre- 
spondence between the absence of nuclear dust and low 
Lfir- For example, of the 19 galaxies in common be- 
tween our sample and that of Lauer et al. (2005), all 
11 of the ellipticals with L70 > 10 40 5 erg s _1 have de- 
tectable nuclear dust, while all but one of the 9 galaxies 
with L70 < 10 40 5 have no nuclear dust detected. Never- 
theless, even for these galaxies with lower L70 we believe 
that most of the Lfir comes from a diffuse component. 
In addition to the argument in Appendix B, our pref- 
erence for extended excess dust is supported by noting 
that all six galaxies with extended FIR emission in Table 
1 (e.g. FW0.1M > 50") have excess FIR emission, i.e. 
L70 > 10 40 - 5 erg s- 1 . 

6.3. External Origin for the Excess Dust 

If the FIR emission in excess of that expected from our 
simple model SED in Figure 11 is from cold dust diffusely 
distributed throughout the hot interstellar gas, where 
has it come from? It is often supposed that optically- 
observed dust in elliptical galaxies has an external ori- 



gin, resulting from a recent merger with a small dust-rich 
galaxy that is otherwise difficult to observe. We adopted 
this explanation previously in our discussion of the ex- 
cess dust in NGC 4636 (Temi, et al. 2003), which we 
interpreted as a dusty merger that occurred within the 
last ~3x 10 8 yrs, the estimated sputtering time in the 
hot gas in this galaxy. However, it has since been de- 
termined that the mean stellar age in NGC 4636 is very 
old (10.3 Gyrs, Sanchez-Blazquez et al. 2006), and this 
strongly argues against a merger ~3x 10 8 yrs ago with 
a dusty, star-forming galaxy. Moreover, we observe that 
the additional dust in NGC 4636 is extended (FW0.1M 
= 54.4" in Table 1) and appears to be distributed rather 
evenly about the center unlike the irregular distribution 
that might be expected to result from a recent merger 
undergoing dynamical friction with accompanying ram- 
stripping and tidal disruption. So our previous interpre- 
tation of NGC 4636 needs to be reevaluated. 

The dust sputtering time for grains of initial radius a in 
the hot gas of an elliptical galaxy is t sp ss a/\da/dt\ sput ~ 
1.2 x 10 8 a AI (n e /10 -2 ) -1 yrs (where is in microns) pro- 
vided the ambient gas temperature is T > 10 6 K. (Drainc 
& Salpeter 1979; Tsai & Mathews 1995). Even if the dust 
is mixed with warm gas at T ~ 10 4 K, where sputtering 
is negligible, its overall sputtering destruction time t sp 
cannot be significantly reduced since, as discussed be- 
low, we generally expect that the warm gas is thermally 
heated and joins the (spatially extended) hot gas phase 
on a time scale ~ 10 6 yrs that is very much less than t sp . 
Therefore it is reasonable to assume that the diffuse FlR- 
emitting dust in our sample galaxies is in direct contact 
with the hot gas during most of its sputtering lifetime. 

The dynamical (freefall) time of a merging dusty 
galaxy across the optical image of a large elliptical is sev- 
eral ~ 10 8 years, comparable to t sp . Therefore, on the 
merger hypothesis, diffusely distributed cold dust must 
be resupplied by mergers in a regular fashion at least 
every ~ 10 8 yrs. Therefore, an entourage of dusty merg- 
ing galaxies should be commonly observed in the outer 
regions of the majority of elliptical galaxies having FIR 
excess emission - but they are not. In fact, mergers be- 
tween elliptical and (small) dust rich galaxies are rarely 
observed at zero redshift. While we believe that mergers 
do occur and that some of the most FIR luminous galax- 
ies in our sample can only be explained by mergers, the 
similarity of the sputtering and dynamical time scales 
makes it difficult to accept the proposition that merg- 
ers are the dominant or only source of (extended) dust 
in most luminous elliptical galaxies. This point of view 
is supported by Figure 5, showing that the incidence of 
young stars in elliptical galaxies is uncorrelated with the 
FIR excess luminosities seen in Figure 3. 

6.4. Internal Origin for Excess Dust 
6.4.1. Evolution of Dust Ejected from Stars 

Before exploring the origin of the additional dust re- 
quired to produce FIR luminosities in excess of L70 <~ 
10 40 ' 5 erg s _1 , it is useful to review the evolution of cold 
dusty gas that is introduced into the hot, X-ray emitting 
interstellar gas by mass-losing red giant stars. Our SED 
model is based on an internal stellar source of dust which 
is confirmed by mid-IR observations of (apparently out- 
flowing) circumstellar dust (Knapp et al. 1989; Athey et 
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al. 2002; Bressan et al. 2006; Bregman et al. 2006). We 
think that the fate of this dusty gas varies with galactic 
radius. 

Throughout most of the volume occupied by galactic 
stars beyond ~ 1 kpc, we assume that both the dust and 
gas ejected from red giant stars ultimately merge with 
the hot gas which is the dominant interstellar phase. 
We expect that gas recently ejected from stars can be 
maintained in pressure equilibrium at warm tempera- 
tures [T « 10 4 K) by photoionization by UV emitted 
by post-AGB stars. Emission from warm gas at r ^> 1 
kpc is faint and has not been observed. Eventually this 
warm gas thermally mixes with nearby hot gas. When 
the dust in this gas confronts the hot gas it will be de- 
stroyed by sputtering in ~ 10 8 yrs. 

However, within ~1 kpc of the centers of elliptical 
galaxies the evolution of dusty stellar ejccta is different. 
Extended warm gas in this region is commonly observed, 
but its velocities differ from those of nearby stars (e.g. 
Caon, Macchetto & Pastoriza 2000). Evidently the hot 
and warm gas exchanged momentum when the warm gas 
had a much lower density. After dusty stellar envelopes 
and winds are ejected from orbiting red giant stars, the 
gas expands until its density approaches (but still ex- 
ceeds) the local hot gas density. At this time the ejected 
gas, which is still moving with the stellar orbital velocity, 
is disrupted by hydrodynamic instabilities and shocked 
to a few million degrees where its density is comparable 
to that of the hot interstellar gas. During this period 
of interaction between the stellar ejecta and the hot gas, 
momentum is exchanged and the ejected gas acquires 
some of the velocity of the ambient hot gas. Even if the 
gas ejected from stars is (transiently) heated completely 
to the temperature of the local hot gas (<~ 10 7 K), its dust 
content will cool the gas rapidly back to ~ 10 4 where it is 
stabilzed by UV from post-AGB stars. This rapid dust- 
assisted cooling occurs because of very frequent inelastic 
collisions between thermal electrons in the hot gas and 
the dust grains. Mathews & Brighcnti (2003) showed 
that efficient dust cooling of the hot gas is also faster 
than the sputtering lifetime of the dust, so the gas retains 
some dust after it cools to ~ 10 4 K. Since the cooling time 
to ~ 10 4 K within 1 kpc is less than the dynamical time 
in the local (stellar) potential, these dense warm clouds 
freefall toward the galactic center and accumulate there. 
As dusty warm gas collects into a disk near the center, 
similar to those observed by van Dokkum & Franx (1995) 
and Lauer et al. (2005), it eventually becomes optically 
thick to the stellar UV radiation. Because of its residual 
dust, we expect molecular formation and eventual cool- 
ing far below 10 4 K. Dense, cold nuclear dust disks and 
clouds similar to those observed can form naturally from 
stellar mass loss in elliptical galaxy cores (Mathews & 
Brighcnti 2003). Dusty galaxy mergers are not required. 
We speculate below that this central dust reservoir may 
be the source of extended diffuse FIR emission in many 
ellipticals. 

Central dusty clouds and disks in ellipticals are not 
merely an amusing curiosity, but are likely to be key to 
understanding the energetics of AGN feedback in reg- 
ulating cooling flows in the hot gas within and around 
elliptical galaxies. Indeed, Simoes Lopes et al. (2006) 
find a correlation between nuclear dust and AGN activity 
in elliptical galaxies. Consider the black hole accretion 



luminosity expected in a large group-centered elliptical 
galaxy by dusty mass loss just within the central kpc 
where rapid dust-assisted cooling is expected. Accretion 
is expected to quickly follow cooling. The AGN luminos- 
ity is L agn = eMc 2 = 6 x 10 45 [Mi/(M Q yr- 1 )](e/0.1) erg 
yr~ 4 where M\ is the total rate of stellar mass loss within 
the central kpc. Since a 10 10 year old stellar population 
of mass ejects 1.5(M*/1O 12 M ) M© yr -1 , the stellar 

mass within 1 kpc, ~5x 10 10 M©, generates Mi w 0.07 
Mq yr _1 . When accreted by the central massive black 
hole, this translates into a luminosity L agn « 4 x 10 44 
ergs s _1 that exceeds the entire X-ray luminosity of the 
most luminous galaxy groups, L x ~ 3 x 10 43 ergs s _1 . 
Consequently dusty stellar mass loss within the central 
kpc of bright ellipticals can generate enough accretion 
energy to balance the entire bolometric X-ray luminos- 
ity of a large galaxy group. The heating of the hot gas 
as a consequence of this accretion also provides enough 
energy to drive buoyant mass circulation in the hot gas 
out to 5-10 kpc. 

The evolution of dusty disks in galactic cores that we 
just described cannot continue very long. If small core 
disks accrete dusty gas at M\ w 0.07 M Q yr _1 , for <~ 10 10 
yrs, and allowing for a larger mass-loss rate in the past, 
currently observed disks should have a typical mass of 
~ 10 10 M© and a dust mass of ~ 10 8 , both far in excess 
of those observed. (Small stellar disks with masses ~ 
10 9 — 10 10 Mq are not commonly observed in the cores 
of most ellipticals, so normal star- formation is ruled out.) 
A large minority of bright ellipticals show no evidence of 
central dusty disks/clouds. Where has this dusty gas 
gone? 

6.4.2. Explosive Ejection of Dust from the Galactic Core 

The excess diffuse dust observed in FIR-emission may 
have been exploded out from the centers of elliptical 
galaxies by the release of AGN energy. We estimate 
from our model SED that the total mass of excess dust 
required to explain the FIR luminosities in Figure 3 is 
M d <~ 3 x 10 5 Mq. This is comparable to the esti- 
mated dust masses in the central dusty disks or lanes 
(e.g. Goudfrooij & de Jong 1995). If the (very uncer- 
tain) dust to gas mass ratio is S « 0.01, then the mass 
of gas that initially contained this dust is M ~ 10 7 — 10 8 
Mq. If cold dusty gas is ballistically ejected from the 
core out to a few kpc, velocities of ~ 100 km s _1 are re- 
quired. The kinetic energy needed to explosively trans- 
port dusty gas from the galactic core to <~ 1 kpc is Ek ~ 
0.5Mv 2 = 10 55 (0.01/£) (M d /10 6 M Q )(i;/100km/s) 2 ergs. 
Although very uncertain, the energy involved is certainly 
larger than that of a single supernova. Since multiple 
starburst activity in elliptical galaxies is very rare, AGN 
energy is required for dust ejection from the core. If the 
AGN dust ejection process is intermittent, this would 
explain those elliptical galaxies that are not observed to 
have optically obscuring dust in their cores or extended 
FIR excess emission. The FIR variance expected from 
intermittent AGN outbursts obviates the argument of 
Forbes (1991) and Simoes Lopes et al. (2006) that the 
failure of an Lfir-Lb correlation necessarily requires an 
external source of dust by mergers. Some explosive char- 
acter to the AGN-dust cloud interaction is implied by 
the chaotic nature of the central dust observed in many 
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HST images (Lauer et al. 2005). This chaotic, denser 
dusty gas is highly transient and must have been ejected 
from the core. These chaotic clouds cannot be under- 
stood as infalling stellar ejecta (as suggested by Lauer 
et al.), formed by the coalescence of many dusty stellar 
envelopes produced by individual stars, since the free fall 
time (~ 10 6 yrs) to the galaxy center is too short for such 
accumulations to occur. 

However, explosive ejections of central dust to ~ 5 — 10 
kpc with energies Ek > 10 55 ergs would produce kpc- 
sized central cavities, strong shocks and heated cores that 
are rarely observed in X-ray images while extended ex- 
cess FIR emission is common. This may be the best 
argument against explosive dust transport. 

6.4.3. Buoyant Ejection of Dust form the Galactic Core 

Density irregularities and multi-temperature X-ray 
spectra of the hot interstellar gas in elliptical galaxies 
reveal buoyant activity, presumably related to a less ex- 
plosive heating of the hot gas with AGN energy (Buotc, 
Lewis, Brighenti & Mathews 2003a, b). In view of these 
X-ray observations, it is plausible that central dust can 
be buoyantly transported out to ~ 5 — 10 kpc, producing 
the extended FIR emission that we observe. 

The buoyant rise time t{, uoy of heated gas is expected 
to be comparable to or somewhat less than the sputter- 
ing destruction time in the buoyant gas. For example 
in the galaxy/group NGC 5044 the freefall time from ra- 
dius rkpc to the galactic center is tff — 4.5 x 10 6 rJ? p 8 c 7 yrs, 
based on a mass model that combines a de Vaucoulcurs 
stellar profile (M* = 3.37 x 10 11 , M Q ; r e = 10 kpc) 
with an Navarro-Frenk- White dark halo (mass M^aio — 
4 x 10 13 M@; concentration c = 8.67). The buoyancy 
rise time tbuoy — Ptff is greater than t ff by some factor 
(5 > 1. ^From Draine & Salpeter (1979) the approxi- 
mate grain sputtering time is t sput ~ 2.4 x 10 6 a /i n e (r) _1 
where is the grain radius in microns. Adopting 
n e w 0.085rj^, cm~ 3 for the observed hot gas density 
in NGC 5044 within 30 kpc, the local sputtering time is 
nearly proportional to tff, t sput « lSa^tff. However, 
the buoyant gas will have a density that is 5 < 1 times 
lower than the observed hot gas density. Therefore, we 
see that the ratio of the local sputtering time in the buoy- 
ant gas to the buoyancy rise time is ~ 1.8a^/ (5/3), which 
may well exceed unity considering the many uncertain- 
ties involved. 

The exact nature of the AGN energy release that initi- 
ates buoyant activity is unclear at present. In about half 
of the ellipticals with optically visible dust, the images 
of the obscuring dust are morphologically "chaotic" and 
distributed in an asymmetric fashion away from the cen- 
ter of the galaxy (Lauer et al. 2005). Evidently the AGN 
energy is sufficient to disrupt the dense central dusty 
disks and lanes which must reform again later by local 
stellar mass loss. It is interesting that the radii of the 
dusty disks observed by Lauer et al. (2005) do not corre- 
late with the stellar rotation, also suggesting short disk 
lifetimes and intermittent disruption of the dusty disks. 

Regardless of the mode of outward dust transport, ex- 
plosive or buoyant, when the dust rises to 5-10 kpc it may 
cause the hot gas to cool back to warm gas temperatures 
~ 10 4 K. The correlations in Figure 8 between F70 and 
i^i6o and the optical line fluxes, if real, indicate a con- 



nection between the warm gas and (excess) dust. Similar 
morphologies of gas and dust images have been discussed 
by optical observers (e.g. Goudfrooij et al. 1994). Al- 
though optical line fluxes appear to correlate with the 
FIR and B-band fluxes, Figure 2 shows that the FIR 
and B-band fluxes do not mutually correlate. This can 
be understood if most of the FIR-emitting dust comes 
not from recent stellar mass loss, but is transported out 
from a reservoir in the galactic core, which contains a 
large mass of dust unrelated to the the global stellar B- 
band flux or luminosity. 

7. CONCLUSIONS 

Our main conclusions are: 

(1) The optical luminosities Lb of elliptical galaxies cor- 
relate with luminosities at 24/im (L24) but not with L70 
or ii6o- This indicates that the transition between cir- 
cumstellar and interstellar dust occurs for dust grains 
with emission that peaks somewhere between 24/z and 
70/rni. 

(2) Elliptical galaxies with similar Lb have L70 and L160 
that vary over a factor ~ 100. This variation is reduced 
to ~ 30 if abnormally dust-rich ellipticals are excluded 
from the sample. 

(3) Neither L^ /Lb nor Li 60 /£_b correlates with the stel- 
lar age or metallicity of the galaxies. 

(4) Weak correlations appear between optical line emis- 
sion fluxes from warm (~ 10 4 K) gas and F70 and .Fi6o- 
This suggests a relationship between dust and gas that 
may have cooled from the hot phase. 

(5) Seven galaxies in our sample are clearly extended at 
70/z with FW0.1M exceeding 50" corresponding to 5-10 
kpc. 

(6) The luminosities and spectral energy distributions 
of elliptical galaxies with the lowest L70 and Liqq are 
consistent with a simple steady state dust model in which 
dust expelled from stars is eroded by sputtering in the 
hot interstellar gas. 

(7) We show in Appendix B that the observed emission at 
70 and 160/xm cannot come from small optically visible 
clouds in the galactic cores similar to those seen in HST 
images. 

(8) In view of the common occurrance of spatially ex- 
tended emission at 70/xm in many ellipticals, it is un- 
likely that this dust can be supplied by successive merg- 
ers with gas-rich galaxies, each occurring during the rela- 
tively short sputtering lifetime of the diffuse dust, ~ 10 s 
yrs. Addiiontal supporting evidence for such mergers is 
lacking. 

(9) The appearance of extended FIR-emitting dust in 
otherwise normal elliptical galaxies can be understood 
by the disruption and buoyant transport of centrally ac- 
cumulated dust by intermittent AGN activity. The buoy- 
ant rise time out to 5-10 kpc, ~ 10 7 yrs, is less than the 
sputtering lifetime of the dust, ~ 10 s yrs. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
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TABLE 1 



Name Sample 3 - T b Flux Density D c R e 70/im FWHM 70/im FW0.1M 





(Pi) 




24 fim 
(mjy) 


70 (im 
(mjy) 


160 fim 
(mjy) 


(Mpc ) 






(Kpc) 


(") 


(Kpc) 


NGC 507 


Temi 


-3.3 




< 2.8 




71.99 


77.1 










NGC 533 


Temi 


-4.8 


7.1±3.1 


23.7±4.9 


26.1±5.6 


68.23 


47.5 










NGC 584 


Kcimicutt 


-4.6 


47.6±7.6 


49.8±6.9 


< 17.2 


23.76 


27.4 










NGC 596 


Suracc 


-4.3 


15.9±4.6 


21.4±5.1 


16.7±7.3 


23.87 


30.0 










NGC 636 


Suracc 


-4.8 


10.3±3.6 


< 16.4 


< 24.5 


22.28 


18.9 










NGC 720 


Temi 


-4.8 


25.6±5.1 


18.2±4.9 


< 21.0 


22.29 


39.5 










NGC 855 


Kcimicutt 


-4.6 


72.3±8.5 


1335±36 


2086±46 


9.33 


18.9 


22.5 


1.0 


58.2 


7.6 


NGC 1377 


Kcimicutt 


-2.1 


1603±40 


4988±71 


2375±48 


22.20 












NGC 1395 


Kancda 


-4.8 


45.4±6.5 


129±11 


194^14 


21.98 


45.4 










NGC 1399 


Temi 


-4.5 


59.9±7.3 


15.9±7.6 


24.1±8.6 


19.40 


42.4 










NGC 1404 


Kcimicutt 


-4.7 


56.6±6.3 


32.6±6.3 


< 32.2 


19.40 


26.7 










NGC 1407 


Kancda 


-4.5 


42.5±6.4 




< 35.4 


22.08 


71.8 










NGC 1700 


Suracc 


-4.7 


18.7±4.9 


29.3±5.3 


40.3±7.2 


38.04 


13.7 










NGC 2325 


Temi 


-4.7 


27.0±6.4 


36.5±6.8 


21.7±6.4 


31.92 


86.5 










NGC 2768 


Temi 


-4.3 


46.6±8.2 


694 ±26 


377 ±29 


22.38 


53.5 


21.4 


2.3 


56.4 


6.1 


NGC 2974 


Kancda 


-4.7 


62.5±12.4 


682 ±18 


1979±54 


30.33 


36.9 


22.8 


3.4 


55.4 


8.1 


NGC 2986 


Temi 


-4.7 


13.7±4.8 


< 2.8 


< 14.5 


32.51 


41.4 










NGC 3265 


Kcimicutt 


-4.8 


287 ±17 


1952±44 


2312±49 


21.33 












NGC 3379 


Fazio 


-4.8 


65.5±8.8 


60.5±8.0 


59.1±8.2 


10.71 


35.2 










NGC 3557 


Temi 


-4.8 


30.5±4.5 


263 ±16 


247 ±18 


34.51 


37.8 


20.0 


3.3 


38.8 


6.5 


NGC 3610 


Suracc 


-4.2 


18.1±4.5 


19.1±7.3 


32.8±16.5 


29.24 


12.8 










NGC 3923 


Temi 


-4.6 


42.9±6.7 


22.7±5.4 


44.0±7.8 


19.14 


53.3 










NGC 3962 


Kancda 


-4.7 


21.5±5.6 


374 ±18 


495 ±24 


23.23 


34.4 


22.6 


2.5 


44.0 


5.0 


NGC 4125 


Kcimicutt 


-4.8 


53.7±6.9 


773 ±28 


1395±37 


27.79 


59.9 


21.6 


2.9 


56.8 


7.7 


NGC 4365 


Cote' 


-4.8 


21.7±4.7 


63.8±8.0 


52.9±7.3 


17.06 


57.2 










NGC 4472 


Fazio 


-4.8 


73.3±8.6 


58.1±7.6 


60.4±7.8 


17.06 


104.0 










NGC 4552 


Kcimicutt 


-4.6 


57.3±7.8 


92.1±10.2 


171 ±16 


17.06 


30.0 










NGC 4589 


Kaneda 


-4.8 


14.9±4.5 


262 ±16 


360 ±21 


26.30 


41.4 


22.6 


2.9 


45.6 


5.8 


NGC 4621 


Cote' 


-4.8 


34.2±6.3 


31.7±5.7 


43.4±6.6 


17.06 


46.5 










NGC 4636 


Fazio 


-4.8 


31.2±5.6 


188 ±12 


169 ±14 


17.06 


101.7 


22.3 


1.8 


54.4 


4.5 


NGC 4649 


Fazio 


-4.6 


106 ±10 


46.3±6.8 


< 27.5 


17.06 


73.6 










NGC 4660 


Cote' 


-4.7 


15.2±4.3 


36.4±6.2 


53.7±8.2 


17.06 


12.8 


20.0 


1.7 


39.0 


3.2 


NGC 4696 


Sparks 


-3.7 


22.9±4.7 


127 ±13 


269 ±21 


39.65 


212.4 


22.6 


4.4 


52.2 


10.1 


NGC 4697 


Suracc 


-4.8 


43.8±6.7 


589 ±25 


755 ±28 


16.22 


75.3 










NGC 4915 


Suracc 


-4.5 


11.5±5.5 


29.2±4.8 


54.2±8.9 


46.98 


7.9 










NGC 5018 


Suracc 


-4.6 


60.5±9.1 


1119±34 


1687±41 


32.36 


24.9 


20.4 


3.2 


39.6 


6.21 


NGC 5044 


Temi 


-4.7 


23.0±5.3 


230±14 


242±18 


32.36 


78.6 


23.1 


3.6 


53.6 


8.4 


NGC 5077 


Fazio 


-4.8 


32.4±5.7 


127 ±11 


141 ±13 


32.36 


25.0 


20.2 


3.2 


49.5 


7.8 


NGC 5322 


Fazio 


-4.8 


41.4±6.4 


455 ±23 


625 ±27 


29.79 


35.2 










NGC 5557 


Suracc 


-4.8 


13.2±5.0 


25.1±10.5 


< 42 


49.02 


28.6 










NGC 5813 


Fazio 


-4.8 


14.9±4.9 


58.4±7.6 


34.8±5.9 


30.78 


48.6 


19.6 


2.9 


49.8 


7.4 


NGC 5846 


Fazio 


-4.7 


38.5±6.2 


102 ±10 


117 ±12 


24.55 


82.6 


20.4 


2.4 


43.1 


5.1 


NGC 5982 


Suracc 


-4.8 


13.5±5.1 


37.2±7.2 


59.9±11.7 


40.18 


29.4 










NGC 6703 


Fazio 


-2.8 


19.9±4.5 


41.4±6.4 


17.9±5.1 


32.06 


23.8 










NGC 7619 


Temi 


-4.7 


9.9±4.7 


< 3.8 


< 14.5 


42.85 


32.1 










IC 3370 


Kancda 


-4.7 


43.8±6.7 


649 ±25 


873 ±30 


45.90 


18.5 


20.5 


4.5 


39.7 


8.8 



a Galaxies are selected from the following Spitzcr observing programs: SINGS Legacy Program, R. Kennicutt (PI); GTO program number 69, 
G. Fazio (PI); GO program number 3403, J. Suracc (PI); GO program number 3506, Sparks (PI); GO program number 3619, H. Kancda (PI); 
GO program number 3649, P. Cote' (PI); GO program number 20171, P. Temi (PI). b The morphological type T is taken from the HyperLeda 
database. Distances arc calculated with Hq — 70 km s _1 Mpc -1 . 
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TABLE 2 



Name 


Log L 24 a 


Log L 70 a 


Log Li 60 a 


Log L B a 


Log L x a 


age 


Z D 




(erg/sec) 


(erg/sec) 


(erg/sec) 


(Lb,q) 


(erg/sec) 


(Gyr) 


(dex) 


NGC 507 




< 40.30 




11.02 


42.76 


8.1±1.9 


261+0 071 


IN V_l V_, JOO 


41.04 


41.19 


40.76 


10.96 


42.29 






IN V.I V_, JOI 


40.95 


40.59 


< 39.66 


10.42 


< 40.15 


2.8+0.3 


U.I I OIU.U4:U 


IN VJ V_y !_> C"J 


40.48 


40.23 


39.65 


10.27 


< 39.66 


2.8+0.4 


n Qfin_|_n not; 


IN V_i V_, UOU 


40.29 


< 40.11 


< 39.82 


10.06 


< 40.17 


4.4+0.6 


n Q7«_|_n noQ 
U.OI u^u.u^o 


NOP 790 

IN V_l V_, 1 ZU 


40.62 


40.10 


< 39.70 


10.44 


40.67 


5.4+2.4 


u.iooiu.uoo 


IN V_l V_, OOO 


40.38 


41.27 


41.00 


8.95 


< 39.83 






NGC 1 377 

IN VJ \y Id 1 ( 


42.48 


42.59 


41.81 


9.60 








IN VJ V_y _1_>J,_/>J 


40.86 


40.94 


40.65 


10.50 


40.95 


7.6±1.4 




NOP 1 W<1 


40.87 


39.92 


39.64 


10.58 


41.69 






NOP 1404. 


40.85 


40.56 


< 39.76 


10.41 


41.25 






IN VJ \y A-'-tVJ 1 


40.83 


< 28.83 


< 39.91 


10.64 


41.06 


7.4±1.8 


n qvq-i-n f)44 


wnn 1 7nn 


40.95 


40.77 


40.44 


10.55 


40.66 


2.6+0.3 


n c;nn+n 049 


NGC 2325 


40.96 


40.71 


40.01 


10.66 


40.76 






NGC 2768 


40.88 


41.68 


40.95 


10.63 


40.44 


10.0±7.0 C 


0.14±0.25 c 


NGC 2974 


41.27 


41.94 


41.94 


10.56 


40.64 






NGC 2986 


40.68 


< 39.61 


< 39.86 


10.57 


41.02 


11.2±1.9 


0.347±0.045 


NGC 3265 


41.69 


42.15 


41.76 


9.46 








NGC 3379 


40.39 


39.98 


39.51 


10.12 


< 39.60 


10. Oil. 1 


0.299±0.036 


NGC 3557 


41.08 


41.64 


41.15 


10.82 


40.64 






NGC 3610 


40.71 


40.36 


40.13 


10.46 


39.89 


1.7±0.1 c 


0.76±0.16 c 


NGC 3923 


40.71 


40.06 


39.88 


10.58 


40.72 


3.3±0.8 


0.623±0.070 


NGC 3962 


40.58 


41.45 


41.10 


10.34 


< 40.28 






NGC 4125 


41.14 


41.92 


41.71 


10.86 


41.0 






NGC 4365 


40.32 


40.41 


39.86 


10.40 


40.31 


7.9+1. 2 d 


0.344±0.065 d 


NGC 4472 


40.85 


40.37 


39.92 


10.96 


41.49 


9.6±1.4 


0.342±0.046 


NGC 4552 


40.74 


40.57 


40.37 


10.35 


40.77 


12.4±1.5 


0.356±0.034 


NGC 4589 


40.53 


41.40 


41.07 


10.39 


40.42 






NGC 4621 


40.52 


40.11 


39.78 


10.38 


40.08 


8.8±l.l d 


0.299±0.050 d 


NGC 4636 

IN V.I V_, lUOU 


40.48 


40.88 


40.37 


10.57 


41.65 


10. 3+1. 3 d 


n 1 Q1 4-0 1 04 d 


NGC 4649 


41.01 


40.27 


< 39.57 


10.79 


41.34 


14 1±1 5 


n QR2+n 029 


NGC 4660 


40.16 


40.17 


39.87 


9.80 


< 39.45 






IN VJ\J t -t\Jcf U 


41.07 


41.44 


41.30 


10.99 


43.23 






NGC 4697 


40.58 


41.33 


40.98 


10.61 


40.18 


8.3±1.4 


n 1 4s+n ri43 


NGC 491 5 


40.92 


40.95 


40.76 


10.34 


< 40.93 






NGC 501 8 


41.32 


42.21 


41.92 


10.63 


< 40.59 






NGC 5044 


40.90 


41.52 


41.08 


10.76 


42.80 






NGC 5077 


41.05 


41.27 


40.85 


10.32 


40.54 






NGC 5322 


41.08 


41.75 


41.42 


10.73 


40.27 






NGC 5557 


41.02 


40.92 


< 40.68 


10.90 








NGC 5813 


40.67 


40.89 


40.20 


10.70 




16.6±2.2 


0.059±0.042 


NGC 5846 


40.88 


40.93 


40.53 


10.72 


41.71 


14.2±2.2 


0.226±0.051 


NGC 5982 


40.86 


40.92 


40.66 


10.59 


41.22 






NGC 6703 


40.83 


40.77 


39.93 


10.43 


< 40.09 


4.8±0.8 


0.354±0.034 


NGC 7619 


40.78 


< 39.99 


< 39.96 


10.64 


41.69 


15.4±1.4 


0.291±0.031 


IC 3370 


41.48 


42.28 


41.94 


10.82 









a Luminosities are calculated with Hq — 70 km s 1 Mpc 1 . b Ages and metallicities are taken from Thomas et al. 2005, unless otherwise stated. c Age 
and metallicity for this galaxy are taken from Howell 2005. d Age and metallicity for this galaxy are taken from Sanchez-Blazquez et al. 2006. 
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Fig. 1. — Grey scale images and contour plots of a representative sub-sample of Ellipticals. Each galaxy is shown in the three MIPS 
passbands centered at 24 (im (left column), 70 fim (center column), and 160 /im (right column). Galaxies are shown in ascending order 
from lowest flux density to highest flux density, including a non detection and sources with extended far-infrared emission. The contour 
levels decrease by factors of 2 from the central peak value. 
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Fig. 1. — - Continued 
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Fig. 2.— 



Comparison of fluxes at 24, 70 and 160/im with optical B-band fluxes. Downward pointing open arrows indicate upper limits. 
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Fig. 3. — Comparison of luminosities at 24, 70 and 160/im with optical B-band luminosities. Downward pointing open arrows indicate 
upper limits. 
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Fig. 4.— 



Comparision of F70/.F160, a measure of dust temperature, with L70. 
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Fig. 5.— 



FIR to optical luminosity (or flux) ratios plotted against apparent stellar age from optical line indices. 
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Fig. 6.— 



FIR to optical luminosity (or flux) ratios plotted against apparent stellar mctallicity from optical line indices. 
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Fig. 7.— 



Comparison of FIR and X-ray bolometric luminosities. 
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Fig. 8. — FIR fluxes plotted against fluxes in optical line emission. ^From left to right the three-panel columns refer respectively to data 
from Goudfrooij (1994a,b), Macchetto et al. (1996), and Sarzi ct al. (2006). Open triangles indicate upper limits in the direction they are 
oriented and open circles are upper limits in both coordinates. 
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Fig. 9. — Comparison of B-band fluxes with optical emission line fluxes in units of 10 14 erg s 1 cm 2 . The optical fluxes are taken 
from (G) Goudfrooij (1994a,b), (M) Macchetto et al. (1996), and (S) Sarzi et al. (2006) Only non-peculiar E and E/S0 galaxies are plotted. 



22 



Temi, Brighenti & Mathews 




Fig. 10. — Surface brightness distriubtions at 70/^m (in arbitrary units) for eight galaxies compared to the point response function 
(PRF) of the MIPS detector. All curves are normalized to 50. at the center and the arrows show the level at 0.1 of the maximum. 
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Fig. 11. — Comparison of 24, 70 and 160^m fluxes with the our model spectral energy distributions. The open triangle is an upper 
limit of -Fi60 for NGC 1404. The solid and dashed lines correspond respectively to emission from grains with original maximum sizes 
a-max = 0.3/im and 1.0/xm. 
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APPENDIX 
NATURE OF FIR-LUMINOUS GALAXIES 

In this appendix we briefly review other observations of galaxies in Figure 3 and Table 2 having the most luminous 
FIR emission. Our objective is to determine if elliptical galaxies with the most luminous FIR emission can be regarded 
as unusual or abnormal and, if so, what is the maximum FIR luminosity of sample galaxies that can be regarded as 
normal. In the following we emphasize observations that are consistent with large FIR emission - evidence of large 
(~ 10 kpc) extended dusty disks, detections of large masses of cold gas or recent star formation, designations as an 
SO galaxy which often contain cold, rotationally supported gaseous disks, strong radio emission, galaxies with nearby 
dusty galaxies or stars that may have contaminated some observations, evidence of a recent merger, etc. - but we 
do not mention attributes that are found in many typical ellipticals: extended optical emission lines perhaps with 
a LINER spectrum, small central (< 1 kpc) dust clouds that obscure starlight, weak central radio sources, X-ray 
emission L x which is highly variable among ellipticals, evidence of old mergers, etc. Galaxies are listed in descending 
order of L 7a . 

IC 3370: Classified as an elliptical in RC2 an RSA, but regarded by Jarvis (1987) and Samurovic & Danziger (2005) 
as an SO or SOpec with an extremely boxy profile. A bright IRAS source in FIR. Richtcr ct al. (1994) detect an HI 
mass of 5 x 10 8 Mq. Contains a prominent dust lane in the inner regions with associated Ha + [Nil] emission and 
possible evidence of ongoing star formation (Michard 2006). Leonardi & Worthcy (2000) and Tantalo & Chiosi (2004) 
suggest an intermediate age population for the central regions. 

NGC 5018: Forms a pair with NGC 5022 7.2' away. Carollo & Danziger (1994) describe a very strong dust lane 
that obscures the innermost regions that may contain recent star formation. They also suggest that the galaxy is very 
flat (about E6). Detected by IRAS at 60 and 100/im. Xilouris et al. (2004) find excess MIR emission at 6.7 and 15/um, 
perhaps suggesting heating by recent star formation. 

NGC 2974: A rapidly rotating E4 galaxy (Emscllcm ct al. 2004). Strong IRAS emission. Kim et al. (1988) 
observes a massive HI disk, Mm = 2 x 10 9 Mq (with D = 38 Mpc). Sarzi et al. (2006) observe a nuclear bar with 
two spiral arms emerging. A bright star is seen in NED images 0.7' to SW. 

NGC 4125: Filho et al. (2002) regard this as an E6pec galaxy. Wickland et al. (1995) detect tentative molecular 
hydrogen emission with mass Mh2 = 7 x 10 7 Mq. Verdoes Kleijn & de Zeeuw (2005) observe a central dust lane of 
size 160 pc (D = 26 Mpc). 

NGC 5322: Carollo et al. (1997) describe a completely obscuring dust lane (perpendicular to the radio jet) that 
hides the center of the galaxy. IRAS detection at 60 and lOO^rn. Proctor & Sansom (2002) and Denicolo et al. (2005) 
find evidence of a ~ 2 — 4 Gyr old young stellar population. 

NGC 2768: Referred to as an E6* in RC3, this is almost certainly an SO galaxy (Funes et al. 2002) with faint 
dust patches surrounding the galactic bulge. Significant dust extinction is visible along the minor axis (Michard 1999). 
Bertola et al. (1992) notes that the ionized gas is dynamically decoupled from the stars, which may be related to gas 
orbiting in a polar ring (Mollenhoff et al. 1992; Fried & Illingworth 1994). Wiklind et al. (1995) observed CO emission 
with an estimated molecular mass of ~ 2 x 10 7 Mq. 

NGC 3557: An E3 galaxy containing an FRI radio source having a double-tailed radio source with a jet and central 
knot (Birkinshae & Davies 1985). Detected by IRAS at 60 and lOO^m. Rampazzo et al. (2005) report a possible dust 
ring near the center and Goudfrooij (1994b) finds an unusual optical emission line image associated with the minor 
axis. 

NGC 5044: Detected by IRAS at 60 and lOO^rn. Very bright and asymmetrically extended optical emission 
(Ferrari, et al. 1999; Goudfrooij et al. 1994a). 

NGC 3962: The Ha+[NII] image shows two sub-systems: an elongated (dusty) central component strongly mis- 
aligned with the stellar isophotes extending to ~ 7 kpc with regular kinematics and a more extended arm-like structure 
associated with bluer stars (Buson et al. 1993; Goudfrooij, et al. 1994b; Zeilinger et al. 1996). Detected by IRAS at 
60 and 100/im. 

NGC 4589: Prominent dust lane near minor axis extending to 20" (Goudfrooij, et al. 1994a, b) Complex gas and 
stellar kinematics suggest a merger remnant (Mollenhoff & Bender 1989). Detected by IRAS at 60 and lOO^m. 

NGC 4697: An E6 galaxy, but may be an SO (Koprolin & Zeilinger 2000). Sofue and Wakamatsu (2003) detect a 
molecular mass of 3 x 10 7 Mq. 

NGC 5077: The E3 designation in RC2 is wrong, this is an SO galaxy (Sandage & Bcdke 1994). 
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ESTIMATED FIR EMISSION FROM CENTRAL OPTICALLY THICK CLOUD 

In this Appendix we determine if a large fraction of the FIR emission from elliptical galaxies can come from small, 
dense, dusty clouds near the centers of the galaxy, or, as we have generally assumed, comes from FIR-emitting dust 
diffusely distributed throughout the inner galaxy. Specifically, can reprocessed starlight from internal and external 
stars heat small dusty regions sufficiently to account for the large FIR luminosities shown in Figure 3? We describe a 
simple calculation to estimate the maximum luminosity at 70^m L70 expected from a small optically thick dust cloud, 
including dust heating by both internal and external starlight. 

Consider a spherical cloud of dusty cold gas with radius r c placed at the center of an elliptical galaxy having a de 
Vaucouleurs luminosity profile. Since the stellar density is greatest at the center, placing the cloud there maximizes 
the energy from internal starlight that can be absorbed and reemittcd at infrared wavelengths. The mean intensity of 
starlight in an (optically thin) elliptical galaxy also is greatest at the galactic center, where the stellar flux incident on 
the dusty cloud from the outside is also maximal. 

Without investigating in detail how the dust is heated and reemitted in the infrared from small, dense clouds, we 
assume that this type of reprocessing can occur. To estimate the fraction of the the FIR that appears in the Spitzer 
70/im pass band, we assume that the overall infrared SED can be approximated with a modified blackbody, with 
emissivity j nu oc v$ B nu (Td) where depends on the properties of the dust and Td is the dust temperature. We 
assume = 1 and Td = (Td) = 27. 9K, the mean values in our SED fits to ISO (Infrared Space Observatory) FIR data 
for elliptical galaxies (Temi et al. 2004). The fraction of FIR emission that emerges in the Spitzer 70fim passband 
(width AA70 = 19/im) is 



L70 _S^ 7 yB u {(T d ))dv 



Lfir / °° V PB v {(T d ))dv 



0.147 (Bl) 



where Aj/7o = cAAyo/A^Q. 

To be specific we consider a typical elliptical in Figure 3 with log is = 10.6. Based on the galaxy-size observations 
of Shen et al. (2003), we assume a typical effective radius for this galaxy, r e = 7.67 kpc. To determine the total 
B-band luminosity of stars inside the dusty cloud, we integrate the de Vaucouleurs profile from the center to radius r c . 
To convert the B-band luminosity to bolometric luminosity of an old single stellar population (SSP), we assume solar 
abundance for the stars, a Kroupa IMF and an age of 10 10 yrs, resulting in (Lb/ Lbq) / (L^oi/ Lboi,o) — 0.396 (private 
communication from Justin Harker). With this information we can estimate the total stellar luminosity produced 
within a dusty cloud of any chosen radius r c and the fraction of this energy that is observed in the Spitzer 70^m 
passband, Lr (int). 

The dusty cloud can also be heated externally by stellar radiation incident on its outer surface. To estimate the 
flux of stellar radiation incident on the dusty cloud we first determine the mean B-band mean intensity within 10-100 
parsecs of the galactic center, Jb ~ 0.03 erg cm~ 2 s _1 ster -1 , found by integrating over the de Vaucouleurs emissivity 
profile corresponding to Lb = 10 10 6 Lbq and r e — 7.67 kpc. The B-band flux on the surface of the cloud is then 
Fb = ttJb- The bolometric flux is larger by a factor F^oi,®/ 'Fb,q = 1/0.396 so the rate that radiative energy from 
external galactic stars is absorbed by the dusty cloud is L-io(ext) — aAirr^iFs / 0.396) where a = 0.5 is the assumed 
albedo. 

Figure Bl shows how the 70/xm luminosities Lyo(int) and L7o(ext) vary with the radius r c of the dusty cloud where 
we see that Lj (ext) < Lj (int). But the important result is that clouds of size r c ~ 100 pc are required to produce 
luminosities L70 ~ 10 41 erg s _1 typical of those plotted in the central panel of Figure 3. In many of the elliptical 
galaxies in Figure 3, a ~ 100 parsec dark cloud would be easily visible even without HST resolution. We conclude 
that it is very unlikely that the large excess FIR emission in Figure 3 above that predicted by our simple model for 
recently produced diffuse dust, can be explained by emission from small dense clouds even under the most favorable 
circumstances. 

More importantly, the approximation of spherical clouds is seriously in error. The gas temperatures associated with 
these dusty clouds (~ 100K?) are much less than the virial temperature of the stars, T» ~ 10 6 K required to support 
a spherical cloud as we hypothesized here. The gas temperature within the dusty cloud is insufficient to provide 
pressure support against the stellar potential due to embedded stars or the external hot gas pressure. Consequently, 
we expect that most of the nuclear dust clouds observed in ellipticals are highly flattened, nearly two-dimensional 
pancakes or disks. In this case we imagine that the L7o(ext) luminosities in Figure Bl are more relevant, further 
reducing the predicted 70/xm emission. We conclude that most of the FIR emission observed at 70fim and 160/z is 
diffusely distributed, not concentrated in dense clouds similar to those seen in absorption by many optical observers. 



26 



Temi, Brighenti & Mathews 




Fig. 12. — Estimated 70/im luminosity emitted from an optically thick dust sphere located a the center of a elliptical galaxy having a de 
Vaucouleurs stellar luminosity with Lb = 10 10 ' 6 Lbq and r e = 7.67 kpc. As the cloud radius r c is varied, the two curves show the 70/im 
luminosity from reprocessed starlight emitted (and absorbed) within the dust cloud Lio(int) (solid line) and the reprocessed luminosity 
absorbed by starlight incident on the sphere from galactic stars beyond r c , L7o(ext) (dashed line). 



